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Abstract —We apply the concept of artificial and controlled 
interference in a triangular relay network with an untrusted 
relay aiming at enhancing the wireless communication secrecy 
between the source and the destination node. In order to shield 
the square quadrature amplitude modulated (QAM) signals 
transmitted from the source node to the relay, the destination 
node designs and transmits artificial noise (AN) symbols to 
jam the relay reception. The objective of our considered AN 
design is to degrade the error probability performance at the 
untrusted relay, for different types of channel state information 
(CSI) at the destination. By considering perfect knowledge 
of the Instantaneous CSI of the source-to-relay and relay-to- 
destlnation links, we first present an analytical expression for 
the symbol error rate (SER) performance at the relay. Based on 
the assumption of an average power constraint at the destination 
node, we then derive the optimal phase and power distribution 
of the AN that maximizes the SER at the relay. Furthermore, we 
obtain the optimal AN design for the case where only statistical 
CSI is available at the destination node. For both cases, our 
study reveals that the Gaussian distribution is generally not 
optimal to generate AN symbols. In particular, the conducted 
numerical results show that the Gaussian distributed AN is 
far from optimal. The presented AN design takes into account 
practical parameters for the communication links, such as QAM 
signaling and maximum likelihood decoding. 

Index Terms —Physical layer secrecy, untrusted relay networks, 
Rayleigh fading, artificial noise jamming, error maximization, 
KKT conditions, square quadrature amplitude modulation, sym¬ 
bol error rate. 

I. Introduction 

T O date, secure data communications heavily rely on 
modern cryptography. Since Diffie and Heilman [2] first 
proposed a key exchange protocol based on computational 
intractability, computational cryptography approaches have 
been extensively studied, such as the popular one based on 
the integer factorization problem [3]. However, the emergence 
of high-performance computers may challenge the existing 
cryptographic algorithms relying on computational hardness. 
As a complementary strategy to provide secure data communi¬ 
cations, physical layer secrecy has recently drawn considerable 
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attention. In particular, physical layer secrecy is viewed as a 
promising solution to provide wireless secrecy in 5G, since 
it does not depend on computational complexity, and has a 
high scalability to allow the coexistence of communication 
terminals with different levels of hierarchical architectures [4]. 
What’s more, physical 1 ayer secrecy can either provide direct 
secure data communication or assist the distribution of cryp¬ 
tographic keys, which makes it particularly favorable in 5G 
networks [4]. Since preliminary works [5,6] characterized the 
secrecy capacity for wiretap channels, secrecy communication 
has been extensively studied for various channel models and 
network setups, such as single-hop wiretap channels [7-12], 
multi-user networks [13,14], and relay networks [15-19]. In 
relay networks, secrecy is an important issue even without the 
presence of external eavesdroppers, since it is often desirable 
for the source to keep its messages confidential from the relay, 
despite the fact that the relay accedes to the request from the 
source to help forward its messages to the destination [20- 
24]. After first proposed in [20], this untrusted relay scenario 
has drawn considerable attention, since it finds diverse and 
important applications in modern communication systems. For 
example, in heterogeneous networks, the relay terminal may 
have a lower secrecy clearance than the source and destination 
pair [22,23]. Therefore, the relay can be partly malicious 
in the sense that it still functions in compliance with the 
relaying protocol, whilst it leaks secret information. Another 
application is the multiple level access control in wireless 
sensor networks, where sensors have different authorizations 
and sensitivities depending on their roles and collected data 
[25-28]. In such setups, the relaying node might be only 
allowed to help forwarding messages from one terminal to 
another, as not all terminals have direct access to other 
members in the sensor group. For instance, in a secure sensor 
network at an airport, terminals with low secrecy levels can 
access limited data, whereas a few terminals at a much higher 
secrecy level are allowed to access all the data [26]. 

To shield the messages from the untrusted relay and enhance 
the secrecy of the wireless communication between the source 
and destination, one popular scheme at the physical layer is 
to introduce controlled artificial noise (AN) to efficiently jam 
the signal reception at the relay. This technique has recently 
been studied from an information-theoretic perspective, such 
as secrecy rate and secrecy outage probability, see e.g., [23, 
29-33]. However, such metrics are in general valid for ideal 
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communication assumptions of continuous input messages and 
random encoding schemes with asymptotically large block 
lengths. In order to take discrete modulation alphabets and 
finite block lengths into consideration, other secrecy perfor¬ 
mance metrics have also been proposed, such as bit error rate 
[34,35]. The observation that in modern wireless communica¬ 
tion systems square quadrature amplitude modulation (QAM) 
is widely used motivates us to address the problem of how to 
optimally apply AN to enhance the physical layer secrecy in 
an untrusted relay network. The symbol error rate (SER) of 
the demodulated signal at the relay is used as a performance 
metric in this paper. 

In this work, we consider the communication between a 
pair of source and destination terminals, where a direct link is 
absent, and therefore, a half-duplex relay terminal is utilized 
to assist the communication [36,37]. To forward the signal 
from the source to the destination, the relay employs the non¬ 
coherent amplify-and-forward (AF) protocol [38], which is 
considered as the most promising solution for current and 
future communication systems, since it offers a reasonable 
tradeoff between the benefits and implementation costs [39]. 
As a result, these relays have already been incorporated in 
Universal Mobile Telecommunications System (UMTS) and 
Release 8 of Long-Term Evolution (LTE) in the form of 
repeaters [36]. For the AF protocol, we make the common 
assumption that the destination has the instantaneous channel 
state information (CSI) of the relay-to-destination link and the 
aggregate source-to-relay-to-destination link [40,41]. Conse¬ 
quently, the destination can obtain the CSI of the relay-to- 
destination link, which is however not available at the relay. 
To secure messages against the relay, we propose a novel 
AN scheme, where the destination designs and broadcasts AN 
symbols to the relay simultaneously with the transmission 
of the symbols from the source. It has been shown that a 
positive secrecy rate can be achievable in this AN assisted 
untrusted relay network [42]. For this setup and with the 
knowledge of instantaneous CSI as assumed in Section III, 
we first investigate the problem of how to optimally design 
the AN to maximize the SER of QAM signals at the relay. 
Note that the requirement of the instantaneous CSI might be 
strong in some practical systems. To reduce the requirement, 
we then investigate the AN design based on statistical CSI in 
Section IV, i.e., the channel variances. For fading channels, 
an important performance metric is the average SER (ASER), 
which quantifies the average decoding error performance over 
fading channels. Interestingly, for both CSI schemes, our 
study shows that it is not optimal to generate AN from a 
Gaussian distribution. For example, based on instantaneous 
CSI, we note that a QAM or a rotated QAM AN generating 
method maximizes the SER at the relay. The results in this 
paper can be used as benchmarks for future analyses of AN- 
based techniques. The main contributions of this work can be 
summarized as follows: 

• By assuming that instantaneous CSI is available at the 
untrusted relay and the destination, upon receiving both 
of the QAM symbols from the source and the AN from 
the destination, an exact expression for the SER at the 


relay when decoding the QAM symbols is first derived. 
Under an average power constraint, this expression is then 
utilized to obtain the phase and power distribution of the 
AN symbols to maximize the SER performance at the 
relay. 

• For the case where only statistical CSI is available, we 
first derive an exact expression for the ASER performance 
at the relay. Next, the optimal power distribution of 
the AN symbols that maximizes the ASER performance 
at the relay is determined. Numerical and simulation 
results demonstrate that the proposed optimal AN designs 
guarantee significant error rate performance enhancement 
compared with conventional AN designs, such as Gaus¬ 
sian distribution. 

Furthermore, the AN-based scheme studied in this paper 
is formulated according to the framework of current cellular 
standards such as LTE/LTE-Advanced of 4G and the next ma¬ 
jor phase of mobile telecommunications standard 5G [4,37]. 
By applying an additional processing unit to generate AN, the 
studied scheme can be easily embedded in a practical system 
to secure wireless communications, e.g., key transmission or 
control signaling. Moreover, our study provides exact SER 
expressions of the QAM signals at the relay, which can be 
used as benchmarks for future extensions, e.g., deriving SER 
expressions for other modulation schemes, and designing AN 
signaling in other scenarios. 

Notation: Throughout this paper, we use CAf{bi, 62 ) to rep¬ 
resent a complex circularly symmetric Gaussian distribution 
with mean bi and variance 62 . We use S {■} and Pr(-) to 
denote statistical expectation and probability. Moreover, | 6 | 
and 6 * represent absolute value and complex conjugate of b, 
respectively. 

II. System Model and Problem Formulation 

Let us consider a relay communication channel as shown in 
Figure 1. A legitimate transmitter (the source) sends informa¬ 
tion symbols to a legitimate receiver (the destination) assisted 
by an untrusted relay. We assume that a direct link between 
the legitimate terminals is not available due to, for example, 
high pathloss. All the terminals are configured with a single 
antenna'. We assume quasi-static Rayleigh fading channels 
and the channel coefficients for the source-to-relay and relay- 
to-destination links are respectively denoted by h and g, 
where h ^ CAf (O, a and g ~ CAf (O, a g). All channels 
are assumed to be reciprocal [43] and constant within the 
transmission duration from the source to the destination. The 
transmission from the source to destination can be partitioned 
into two time slots. During the first time slot, the source 
broadcasts the symbol m. Meanwhile, to prevent the relay 
from deciphering the information, the destination sends the 
AN symbols z to increase the noise level at the relay. For 
this purpose, all the terminals are assumed to be perfectly 
synchronized. Similar to physical-layer network coding [44], 

'Note that this work can be generalized to the scenario where at least one 
of the source and destination has multiple antennas. Whilst beyond the scope 
of this work, the corresponding precoding design at each terminal remains an 
interesting topic for future extension. 
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Fig. 1. The considered system model in this work. 


the synchronization is an important issue, which requires 
further investigation. At the end of the first time slot, the 
received baseband signal at the relay can be mathematically 
expressed as 

y = hm + gz + n (1) 


where the additive white Gaussian noise (AWGN) symbol n 
is assumed to be zero mean and have the one-sided power 
spectral density Nq in Watts/Hz. During the second time 
slot, the source is muted, and the relay applies the amplify- 
and-forward protocol [38] to forward a scaled version of its 
received signal to the destination. Upon receiving this signal, 
the destination removes the AN z prior to decoding m. For 
the afore-mentioned communication protocol, the destination 
is required to have the perfect CSI of h and g in order to 
maximize the system’s transmission rate [45]. To be more 
specific, the destination performs the channel estimation of the 
relay-to-destination link and the compound source-to-relay-to- 
destination link, and the channel coefficient of the source-to- 
relay link is then obtained [41]. 

While exploiting the untrusted relay to help the data trans¬ 
mission, the system is designed to ensure that the untrusted 
relay cannot decipher the source’s symbol m. In addition, 
we assume that the relay can obtain knowledge about the 
modulation scheme of the source, e.g., by tracking the com¬ 
mon control channel of the network. Moreover, it is important 
to note that the AN z in Eq. (1) is solely designed by the 
destination so that the relay is not aware of the existence of 
the AN, neither the AN distribution. Therefore, when decoding 
the symbols transmitted by the source, the relay treats the AN 
from the destination as noise. Furthermore, we assume that 
the source transmits demodulation reference signals (DM-RS) 
[37] so that the relay can perfectly estimate h to perform ideal 
coherent demodulation as 


y = j^y = rn- 


h*g h* 


( 2 ) 


Forwarded by the relay, the destination also has access to h 
and the DM-RS of the source in order to demodulate the signal 
TO. For the source’s symbol to, we consider square M-QAM 
modulation types with M = 4^ and k = 1,2,..., which are 
frequently used in the current and upcoming communication 
standards [37]. We also make the common assumption that 
constellation points are uniformly distributed [46]. In addition, 
we use Em and to denote the average energy per symbol 
for the source’s signal to and the AN z sent by the destination, 
respectively. Finally, Em is assumed to be known at the 
destination for the AN design. 


Since the destination has perfect knowledge of the instan¬ 
taneous channel coefficients h and g, this knowledge can be 
efficiently used to design the AN symbols z. The relay uses the 
ideal coherent demodulation with minimum distance detection 
to recover the source’s signal to, and the corresponding SER 
is the performance metric used in this paper. The objective 
is to find out the AN symbols z at the destination in order 
to maximize the SER performance of the square M-QAM 
modulation at the relay, which we address in Section III. To 
reduce the complexity of the AN design, the AN symbols z 
can also be obtained based on long-term CSI of h and g, i.e., 
and CTg. In this case, the AN symbols z can be determined 
to maximize the ASER at the relay, and this is investigated in 
Section IV. 

III. AN Maximizing Relay’s SER Performance 

The SER analysis in this section consists of three parts. 
In Section III-A, for given instantaneous channel realizations 
h and g, we derive the SER expression at the relay for a 
given AN symbols z of square M-QAM modulation. Based on 
this expression, in Section III-B, we further study the problem 
of how to select the phase of z given its amplitude |z|. The 
problem of assigning power to the AN symbols z is considered 
in Section III-C, where we derive the optimal distribution of 
the power of z, i.e., |zp, based on an average power constraint 
E,. 

A. SER Expression for a Given z 

The SER performance of the square M-QAM modulated 
signaling over AWGN [46] and over fading channels [47] 
has been widely studied. To provide a comprehensive study 
on the SER performance at the relay, we first investigate the 
AN design using the instantaneous CSI in this section, and 
our SER analysis is based on the procedure presented in [46, 
Chapter 5]. 

As an example, the constellation diagram for the square 
16-QAM modulation is illustrated in Figure 2. Denote the 
minimum distance between two constellation points as 2a 
(a > 0), then for a general M-QAM constellation, the average 
energy per signal symbol Em can be expressed as [46, Eq. (5- 
2-76)] 

Em = ^a^Tm (M - 1) (3) 

where Tm denotes the symbol duration. 

Since the AN design in this section is based on h and g, for 
simplicity of notation, we define s := and n := in 
Eq. (2), and denote the real and imaginary parts of s and 
h as Sr, Si and hr, hi, respectively. Thus, we have that 
s = Sr + jsi and h = hr + jhi with j = \/—l. Note 
that the equivalent noise h is identically distributed with n, 
and thus h ^ CMf), 2tT^), where cr^ = ^ vertex 

such as the point “0” in Figure 2, the received symbol y 
in Eq. (2) lies outside the decision region of the point “0” 
if Sr + hr > \h\a or Si + hi < —\h\a. Thus, given s, the 
error probability of the point “0” is given by Pe,o = 1 ~ 
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Fig. 2. The constellation diagram of square 16-QAM with Gray mapping. 


(1 — Pr (sr +nr > \h\a)) (1 — Pr (s^ + fii < —|/i|a)), which 
can be computed as 


Pe,0 = l- 1-Q 


\h\a—Sr 


1-Q 


|/i|a+Si 


(4) 


where Q {■) is the Gaussian Q-function [47, (4.1)]. Sim¬ 
ilarly, the error probability for the ith constellation point 
Pe,i (I = 0,1,..., M — 1) is obtained by computing the 
probability that y in Eq. (2) lies outside its decision region. 
Based on the assumption that the constellation points are 
uniformly distributed, for a given z, or equivalently for a 
given s, the SER of the square M-QAM signal at the relay 
can be derived by averaging the respective error probability 
expressions for all symbols under the assumption of a uniform 
symbol distribution, resulting in 


SER(s) =c 
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Vm-1 


. Using Eq. (3), we have that 




G Y iVo (M - 1) ■ 


(5) 


( 6 ) 


Note that in the case of z = 0, the SER expression in Eq. (5) 
coincides with that of the AWGN channel [46, Eq. (5-2-79)]. 

Before proceeding, it is interesting to note that the SER 
expression in Eq. (5) depends on the real and imaginary parts 
of s, which motivates us to express s as 


s = |g||z|exp(j6») 


(7) 


where we have expressed z, g and h as z = IzlexpQB^), 
g = IpI exp(j0g), h = \h\exp{j9h), respectively, and 6 = 


0g — dh + (Iz- Now, inserting Eq. (7) into Eq. (5), the SER 
expression can be rewritten as a function of 9 and |z| as 


SER(6», |z|) = 
c 


<31 (jfiik+Mlncose 
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• ( 8 ) 


Erom Eq. (8), we can observe that the channel gains \h\ and 
Ip I play an important role on the SER performance, which we 
summarize in the following propositions. 

Proposition 1 (SER decreases in \h\}: The SER in Eq. (8) 
is a monotonically decreasing function of \h\. 

Proposition 2 (SER increases in \g\): The SER in Eq. (8) 
is a monotonically increasing function of |p|. 

To prove these propositions, we denote 


?(|/i|,l5l)=EQ(—+ (-l)'®cos0) (9) 

1=0 ^ ^ ^ ^ 

(10) 


We first focus on Proposition 1. Since |p|) is a de¬ 

creasing function of \h\, the first derivative of 11?!) with 
respect to \h\ satisfies < 0, and for a given |p|, 

Cd^l, Isl) < ■f(0, IffI) = 1- Similarly, < 0 and for a 

given |p|, 77(|ft,|,|p|) < 1. In addition, since 0 < c < 1, the 
first derivative of the SER expression of Eq. (8) with respect 
to \h\ can be written as 


c{l - c^{\h\,\g\)) 
< 0 


M\h\, |g|) 

d\h\ 


+c(l-c77(|/i|,|p|)) 


d\h\ 

( 11 ) 


which proves Proposition 1. This proposition agrees with 
the intuition that a stronger source-to-relay link improves the 
decoding performance at the relay. 

Regarding Proposition 2, for a given \h\, we use the first 
derivative of the Gaussian Q-function 


dQ{x) 1 


dec 


-tt 


( 12 ) 


to obtain 


d^d/i|,l5l) ^cos0 


d|5l 




exp 
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(13) 
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which indicates that |(/|) is an increasing function of |(/|. 
Thus, we have 

C(|/i|,l5l) < lim C(|/j|,l5l) = 1- (14) 

|g|->-oo 

Therefore, since 0 < c < 1, the first derivative of the SER 
expression of Eq. (8) with respect to |p| can be derived as 

f-nu] I l^^ l^?l) I ^-1 /\i^\ I 1'^'^ IpI) 

c(1 - c^{\h\, |g|))- — -+c(l - cri{\h\, IpD)- — - 

> 0 (15) 

which proves Proposition 2. This proposition suggests that a 
stronger relay-to-destination link deteriorates the SER perfor¬ 
mance at the relay. 

A careful observation of Eq. (8) shows that the AN design 
can be decoupled into the optimal design of the phase and 
the amplitude, respectively. In the following subsection, we 
first derive the optimal rotation angle 9 that maximizes the 
SER given by Eq. (8). The value of 6^, which denotes the 
corresponding phase of the AN defined following Eq. (7), is 
then determined accordingly. 


B. Selecting 9z for a Fixed \z\ 

The focus of this subsection is to determine the optimal 
phase distribution of the AN symbols z. In this case, the power 
|zp of the AN, is assumed to be fixed. Then, for simplicity, 
we can omit \z\ in the argument of SER(0, |z|) in Eq. (8) and 
denote the SER expression as SER(0) in the remainder of this 
subsection. 

Due to the ^-periodicity of the function SER (9) and due 
to its symmetry, we only consider the interval 9 G [O, 
to determine the maximum of SER(0). We first numerically 
verify that the following result. 

Result 1: SER(0) is a quasi-convex function of 9 for all 
9 G [O, j] , and either 9 = 0 or 9 = j maximizes the value 
of SER(6»). 

To illustrate this result, consider the constellation point “5” 
in Eigure 2 as an example. Given the power of the AN 
symbols z, in order to maximize the expected SER with 
square 16-QAM modulation at the relay according to Result 
I, the destination can either allocate all the power \z\^ to the 
direction of “12” or “13”, or equally distribute this power 
between the directions of “12” and “13”. Due to the similarity 
of this AN phase design to the constellation diagram, we refer 
to the case of 0 = 0 as the rotated QAM phase selection 
and the latter case of 6* = ^ as the QAM phase selection. 
Accordingly, the optimal phase selection of the AN z is 
computed as 9z = 9h - 9g or 9z = j + 9h - 9g. 

Next we address the question in which case using the rotated 
QAM constellation for the AN maximizes the SER at the 
relay and in which case the M-QAM constellation is optimal. 
Towards this aim, we obtain Eigure 3, which shows |z|/(t as 
a function of a/ a. Each curve in this figure is plotted by 
numerically solving the equation SER(7r/4) = SER(O) for 
|z|/cr for given values of a/ a. The QAM region and the rotated 
QAM region denote the regions of SER(7r/4) > SER(O) 
and SER(7r/4) < SER(O), respectively. Therefore, Eigure 3 


displays phase selection thresholds for different system pa¬ 
rameters. In this figure, a/a and |z|/ct indicate the energy of 
the source signal m and the energy of the AN z, respectively. 
We focus on the relation between the signal energy and the 
AN energy, and thus ignore the effect of channel gains by 
assuming h = 1 and g = 1- 

The curves in Eigure 3 show that for a given signal power, 
there exists a threshold for the AN power, under which the 
optimal phase selection is 0 = 0, and above which the optimal 
phase selection is 9 = tt/A. To explain the intuition behind 
this, consider the constellation point “5” in Eigure 2. With the 
purpose of maximizing the SER performance at the relay, if the 
power level of the AN |z|/ct is low (as compared to a given 
signal power Em), the AN transmitted from the destination 
may deviate the equalized information signals received at 
the relay towards the adjacent constellation points “1”, “4”, 
“13”, and “7”, whilst if the power level of the AN |z|/cr is 
sufficiently high, the AN may shift the receive symbol towards 
the farther constellation points “0”, “12”, “3”, and “15”. Erom 
Eigure 3, one can also observe that for small values of a/a, 
the values of the threshold are very close to zero. This fact 
indicates that QAM phase selection is the preferable scheme 
if the received signal-to-noise ratio (SNR) at the relay is low. 

In particular, from Eq. (8), we have 


lim SER (0, |z|) =c-f 2c(l — c)(5 —L 

l^l^oo y (7 

s/M-l 2{y/M-l) 1 

- M ^2 

M - 1 


M 


(16) 


and 


lim SER (7r/4, |z|) = 2c - ^ ^ (17) 

|z|—>-00 IVl 


where the equality in (16) holds only when a = 0, i.e., when 
no signal is transmitted. This indicates that when the signal is 
present, given sufficiently high AN noise power, the QAM 
phase selection yields a higher SER. This analysis can be 
confirmed by Eigure 3. Also note that is the SER when 
the relay does not have any prior information and randomly 
guesses the value of signal m for decoding, which serves as 
an upper bound on the SER of any AN scheme. In this work, 
we term this as non-informative error performance. Moreover, 
Eqs. (16) and (17) reveal that when the AN power is suf¬ 
ficiently high, theM-QAM selection asymptotically achieves 
the non-informative error performance. 


C. Assigning \z\ for a Given Ez 

In Subsection III-B, we have shown how to select the phase 
9z of z if the amplitude |z| is a fixed value. In general, 
however, it may be optimal to assign different powers to 
different AN symbols for a given average AN power. In this 
subsection, we assume that the average energy per AN symbol 
is limited to an average symbol energy E, i.e., 

Ez = Tm£{\z\^} <E. 


(18) 
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Fig. 4. SER performance of 16-QAM modulation as a function of l-^l/cr, 
where a/cr = 4, h = 1, and g = 1. 


Furthermore, we denote P = E/Tm as the average power of 
z. Based on the results in Subsection III-B, the expected SER 
obtained from optimal phase selection is given by 

^(|z|) = max{SER(7r/4, |z|),SER(0, |z|)} (19) 

where the SER(0, jzl) is defined according to Eq. (8). Eol- 
lowing the proof of Proposition 2, one can show that the 
first derivative of SER(0, |z|) in Eq. (8) with respect to 
\z\ is non-negative, and thus SER(0, |z|) is a monotonically 
increasing function of jz] for a given 9. Accordingly, SER (jzl) 
in Eq. (19) is also a monotonically increasing function of \z\. 
Eigure 4 displays the values of SER(|z|) for 16-QAM and 
a/a = 4, where both the analytical and the simulated results 
are plotted. From the figure, we can see that as |z|/cr increases, 
the rotated QAM and the QAM phase selections alternatively 
achieve a better SER performance. In this figure, note that 
|zp = P can be viewed as a deterministic power usage. 

The objective of this subsection is to derive the optimal 


distribution of |z|^ to maximize the expected value of the SER 
in Eq. (19), given the instantaneous CSI h and g. Alternatively, 
the power needs to be smartly allocated to yield an upper 
bound on the SER performance with deterministic power 
usage. Denoting the probability density function (PDF) of \z\^ 
as /(•), by taking into account the power constraint in Eq. (18), 
the following optimization problem can be formulated 

pCO 

max / SER (v/ir) / (a;) da; (20a) 

fi^) Jo 

pCO 

s.t. / xf {x) dx < P (average power) (20b) 

Jo 

f{x) da: = 1 (total probability) (20c) 

/ (a:) >0, for a: > 0 (non-negativity) (20d) 



where the average power constraint (20b) follows from (18). 
For computational tractability, we only consider PDFs f{x) 
for which all integrals in (20) exist. 

The following theorem provides interesting insights into the 
power allocation problem and meanwhile can greatly simplify 
the computation of the PDF of \z\. 

Theorem I (SER-maximizing power distribution): Any PDF 
f{x) solving the problem (20) has the form 


f{x) = {l-p)5{x-xi)+p6{x-x2) (21) 


where <5 (a:) is the Dirac delta function defined by [48, Chapter 

2 ] 


- T)g{t)dt = g{T) 


( 22 ) 


for any continuous function g and any value of r, and 


P = 


P — Xi 

X2 — Xi 


(23) 


with 0 < a;i < P < a; 2 . 

Proof: See Appendix. 

This theorem means that only two kinds of AN symbols are 
generated from the destination, one with probability (1 — p) 
and one with probability p, and the corresponding powers are 
xi and X 2 , respectively. 

Inserting Eq. (21) into Eq. (20), the expected SER in 
Eq. (20a) can be computed as 


SER(a:i, X 2 ) 

pOO 

= / SER {sfx) f (x) dx 
Jo 

= (1 — p)SER (v^) + pSER {y/xf) 

= ^ ~ sliR {^). (24) 

X2 - Xi X2 - Xi 

Erom Eq. (24), at high average AN power, i.e., P ^ 00 , we 
observe that the expected SER in Eq. (24) approaches 
which corresponds to the non-informative error performance. 
The proof can be sketched as follows. From Eqs. (16) and (17), 
we know that lim| 2 |_^oo SER (|z|) = . Given P —>■ 00 , 

by applying the Optimal Power Allocation Theorem, we have 
X 2 —>■ 00 due to X 2 > P. In this case, if Xi —>■ 00 , we have 
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linia^i^oo SER (y/jrT ) = and lim 2 ; 2 _>oo SER {y^) = 
thereby SER(xi,a; 2 ) in Eq. (24) approaches If 

xi is a finite value, X 2 ~ P, then the term SER (v^) 
in Eq. (21) approaches zero, and the term J^^^SER (v^) 
in Eq. (21) converges to , which explains the above 
observation. 

Another interesting observation is that the maximum 
SER(a:i,a; 2 ) is a monotonically decreasing function in \h\. 
This is because SER(a;i,X 2 ) in Eq. (24) is a linear combi¬ 
nation of SER (y^) and SER(y/T 2 ), each of which is a 
monotonically decreasing function of |/i| as previously shown. 
Eurthermore, maximizing SER(a;i, a; 2 ) with respect to xi and 
X 2 preserves the monotonicity. Similarly, it can be shown that 
SER(a:i, a; 2 ) in (24) is a monotonically increasing function of 
|p|. Therefore, by applying the optimal AN design, the SER 
at the relay increases with the relay-to-destination link quality 
and decreases with the source-to-relay link quality. 

The optimal values of xi and X 2 maximizing SER(a;i,X 2 ) 
in Eq. (24), denoted by x\ and X 2 , can be computed numer¬ 
ically based on h, g, P, M, a, and a. Substituting x\ and 
X*2 back into (24) yields the maximum expected SER, which 
is denoted as SER^ax- For example, in the case of ft, = 1, 
5 = 1, P = 3.9811, M = 16, a = VlO, and cr = 1/^2, the 
optimal values can be computed as xi = 0 and X 2 = 13.7098, 
yielding SER (yil) = 0 and SER {^) = 0.5832. The 
corresponding maximum expected SER is SERmax = 0.1694. 
Interestingly, if we use a deterministic power usage P, i.e., the 
AN PDE can be written as f{x) = 5{x—P), the corresponding 
expected SER is 0.0371. Therefore, applying the Theorem I at 
the destination yields a SER increment of 357% at the relay 
compared with the deterministic power usage. 

Here, we summarize the procedure to design and generate 
AN symbols as follows: i) compute the amplitudes of the AN 
symbols by numerically maximizing Eq. (24), ii) determine the 
corresponding phase value for each amplitude using Eq. (19), 
iii) generate the AN symbols by the obtained amplitudes and 
phases. Erom the procedure, we can see that the computational 
complexity of the AN generation mainly stems from the opti¬ 
mal amplitude computation in step i). There exists numerous 
non-linear optimization algorithms to solve this optimization 
problem, such as the Newton-Raphson method, Nelder-Mead 
method [49], which is known as “fminsearch” implemented in 
MATLAB. Eor example, using “fminsearch” with the default 
termination tolerance, i.e., 10 “^, the iteration number to solve 
the problem is generally under 100 for the numerical examples 
in Section V. Moreover, note that the computation in steps i) 
and ii) is required after the change of the channel coefficients 
ft or g. 

D. Relation Between Phase Selection and Power Allocation 

Having determined the optimal phase selection and the 
SER-maximizing power distribution, it is interesting to exam¬ 
ine their interplay in the SER performance, which we illustrate 
in Eigures 5-7. Eigure 5 depicts the SER performance using a 
deterministic power level |zp = P and the optimal power al¬ 
location. The line denoted as “Non-informative” represents the 
non-informative SER performance. In the case of 16-QAM, the 


non-informative SER is equal to 15/16. Comparing Eigure 4 
and Eigure 5, we observe that for a given SNR E^/Nq, the 
derived Theorem I yields an upper bound on the SER with the 
deterministic power usage, i.e., the curves SER ^0, and 

SER ^ j, s/pj. To take a deeper look, we depict in Eigure 6 
the relation between the optimal power allocation x\ and 
The corresponding phase selection and are obtained 
from evaluating SER(iy7^) and SER(y^7f) in Eq. (19), 
respectively. Eigure 7 plots the probability to transmit i.e., 
p given in Eq. (23), noting that the probability to transmit x\ 
is given by (1 — p). 

In correspondence to Eq. (21), we observe from Eigures 6 
and 7 that as the power of the AN P increases, either the 
probability p remains constant and the powers x\ or, 
increase or the values of x\ and a;^ remain constant (with 
X 2 > x\) and the probability p of transmission with the larger 
power X 2 increases. To simplify the discussion of Eigures 5-7, 
we partition the values of the AN-to-natural noise ratio (ANR) 
per symbol into multiple regions by introducing the transition 
points A-D. Erom these figures, we obtain the following 
observations. 

Region / (0 < E^/Nq < A): In this region we observe that 
when the ANR is low, the optimal AN design is to transmit a 
noise symbol with power a;^ and rotated QAM phase 6 = 0 
at a constant probability po, and to transmit no AN (a;* = 0) 
with probability (1 — po). Thus, as P increases, the power 
X 2 increases while the probability of transmitting AN remains 
constant. 

Region II (A < E^/No < B): When the ANR is medium, 
the power to transmit rotated QAM AN symbols a;^ reaches 
a constant, and the probability to transmit the AN symbols 
increases as E^/Nq increases. The corresponding SER in 
Eigure 5 is higher than those using the deterministic power 
P. 

Region III (B < Ez/Nq < C): In this region, the 
probability to transmit AN symbols with the rotated QAM 
phase selection reaches 1 , and the corresponding transmitting 
power reaches the highest a;^ = P- In other words, for the AN 
symbols, the deterministic power usage and the rotated QAM 
phase selection achieve the maximum SER, which can also be 
observed in Eigure 5. 

Region IV (C < Ez/Nq < D): The ANR arrives at a 
threshold Ez/Nq = C, above which the optimal AN design is 
to use both of the two kinds of AN symbols: the rotated QAM 
phase selection with power x\ and the QAM phase selection 
with power a:^- The probability to transmit the QAM phase 
selection symbols increases as the ANR increases. 

Region V (Ez/Nq > D): When the ANR is large, the 
QAM phase selection with a deterministic power P achieves 
the maximum SER. In particular, when Ez/Nq is sufficiently 
large, i.e., Ez/Nq > 17 dB in Eigure 5, the SER converges to 
the non-informative SER, which was mathematically proved 
in Eq. (17). 

To intuitively understand the above observations, again take 
the constellation point “5” in Eigure 2 as an example. When 
the AN power is low, the optimal AN generation scheme is 
to burst the limited power to move “5” towards the adjacent 
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Fig. 5. SER performance of 16-QAM modulation as a function of Ez/Nq 
without and with optimal \z\, where Em/N q = 20 dB, h = 1, and g = 1. 



Fig. 6. Optimal amplitude distribution and the corresponding phase selection 
for 16-QAM modulation, where Em/N^ = 20 dB, h = 1, and g = 1. 


points “1”, “4”, “13” and “7” in order to introduce decoding 
errors. As the AN power goes sufficiently large, some of the 
power can be used to move “5” towards the points “0”, “12”, 
“3” and “15” in order to induce more decoding errors. Until 
when the AN power is significantly large, beaming all the 
power to the directions of points “0”, “12”, “3” and “15” yields 
the maximum SER. 

IV. AN Maximizing Relay’s ASER Performance 

In Section III, we have obtained the optimal phase and am¬ 
plitude distribution of the AN based on the instantaneous chan¬ 
nel knowledge of the source-to-relay and relay-to-destination 
links. In this section, we consider the optimal distribution 
of the AN symbol powerjzp that maximizes the ASER per¬ 
formance of the M-QAM modulation at the relay provided 
that the long-term statistical CSI is known to the destination. 
Eollowing a similar approach to Section III, we first present 
an analytical expression for the aforementioned performance 
metric over Rayleigh fading channels. The optimal phase and 



Fig. 7. Probability to transmit for 16-QAM modulation, Em/No = 20 
dB, h = 1, and g = 1. 

power distributions of the AN that maximizes this performance 
are then determined. 

A. ASER Expression for a Given z 

In this section, the channel coefficients h and g are assumed 
to be random variables. Moreover, the envelope of the relay- 
to-destination channel coefficient g is assumed to be Rayleigh 
distributed. Eor a given AN symbols z, the received AN 
at the relay gz can be easily shown to be distributed as 
gz ^ CA/" (O, (Tg|zp). This indicates that the received AN 
at the relay is an extra source of the AWGN n. Hence, as 
can be observed from the signal model in Eq. (1), the overall 
received noise at the relay is a superposition of the AN and the 
natural AWGN. The relay-to-destination channel coefficient g 
and the AWGN symbols n are independent, and it can be 
easily shown that the received noise gz + n is distributed as 
gz + n CM (O, CTglzp-f 2cr^), where 2(t^ is the variance 
of the AWGN n. Eor Rayleigh faded source-to-relay channel 
gain I h \ and the overall noise symbol gz + n, the performance 
of the M-QAM modulation at the relay can be obtained using 
[47, Eq. (8.106)] as 



where ^(z) = ^.2 |^'| 2 ^o .2 represents the average received 
SNR per symbol. 

Similar to Propositions 1 and 2 made on Eq. (8), the effects 
of channel statistics on the ASER given by Eq. (25) can be 
summarized as follows. 

























































9 


Proposition 3 (ASER monotonicity): The ASER in Eq. (25) 
is a monotonically increasing function of a^, as well as a 
monotonically decreasing function of 

This proposition is intuitive and is similar to Propositions 
1 and 2 in Section III: Stronger relay-to-destination channels 
help to deteriorate the ASER performance at the relay, whereas 
stronger source-to-relay channels improve the average decod¬ 
ing performance at the relay. 

Proposition 4 (high signal and AN power performance): At 
high SNR and high ANR, i.e., Em ^ cr^ and \z\^ ^ 
is{z) « 

Proposition 4 reveals that provided an adequately high 
signal and AN power is available, the relative strength of the 
source-to-relay and relay-to-destination links affects the ASER 
performance of the QAM signals at the relay. 

Proposition 5 (high AN power performance): In the high 


AN power regime, i.e., for \zf —>■ oo, 7 s —>■ 0, we have 
limi^i^oo ASER(2;) = 2c - 

This proposition shows that given sufficiently high AN 
power, the ASER at the relay is close to the non-informative 
error performance. 

Eurthermore, we observe that the ASER expression in 
Eq. (25) is independent of the phase of the AN, which is 
different from the AN design in Section III in the case of 
instantaneous CSI. Eollowing a similar approach as in Section 
IITC, our objective in the following subsection is to derive 
the SER-maximizing power distribution for the AN under an 
average power constraint. 


B. Assigning \z\ for a Given Ez 

We rewrite the function ASER( 2 ) in Eq. (25) as 
ASER(y^) to reflect the effect of AN power, where x = | 2 :p. 
Similar to Subsection III-C, the optimal power assignment 
problem is given by replacing SER(y^) in Eq. (20) by 
ASER(y^), and now f{x) represents the PDE of x for 
Rayleigh fading channels. Since one can see that ASER(y^) 
in Eq. (25) is a monotonically increasing function of x, the 
distribution of the optimal AN directly follows from the PDE 
expression derived in Eq. (21). Therefore, substituting Eq. (21) 
into Eq. (25) yields the expected ASER 


ASER(a;i, X 2 ) 

= ASER(ViT) + ^ ~ ASER(v^) (26) 

X2 — Xi X2— Xl 

where P = -^ is the average power constraint on the 
AN, and 0 < xi < P < X 2 - The optimal values of 
Xl and X 2 maximizing ASER(a:i, 0 : 2 ) in Eq. (26) are then 
derived numerically based on P, M, Em, cr, and cr^, 
where and cr^ represent the statistical/long-term CSI. Note 
that the computational complexity of the derivation has been 
investigated in Section III.C. 

Similar to the analysis following Eq. (24), we can also show 
that as the AN power increases, the expected ASER in Eq. (26) 
converges to the non-informative SER performance. 
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Fig. 8. SER performance of 4-QAM modulation as a function of E^INq, 
where Em/No = 10 dB, h = 1, and 3=1. 


V. Numerical and Simulation Results 

In this section, we present simulation results to validate the 
obtained analytical results derived in the previous sections. In 
all figures, the average SNR is chosen as Em/Ng = 10 dB. 

In order to depict the SER performance at the relay using the 
optimal AN design in Section III, we plot Eigure 8, where the 
signal symbol m is randomly selected from the square 4-QAM 
constellation. We observe that all analytical results coincide 
very well with the corresponding simulation results. Eor the 
sake of comparison, the expected SER values for Gaussian 
AN are also plotted in Eigure 8. In the Gaussian case, we 
assume that the AN symbols 2 is generated according to a 
complex Gaussian distribution with the same average symbol 
energy Ez- The line denoted as “Non-informative” represents 
the non-informative SER performance. The curve “w/o AN” in 
the figure depicts the SERs at the relay without AN, which lies 
a value of 0.0016 for the given parameters. Eigure 8 shows 
that the SER at the relay can be significantly increased by 
applying our AN design. Even when the ANR Ez/Nq is small, 
e.g., Ez/Nq = 2 dB, the SER is increased from 0.0016 to 
0.05, i.e., an increment of 303%. Moreover, Eigure 8 clearly 
demonstrates that the Gaussian distribution is not optimal 
for our AN symbol generation and our scheme described in 
Section III can induce much larger SER at the relay. Eor the 
case shown in Eigure 8, if the ANR Ez/Nq is above 16 dB, the 
maximum SER, i.e., SERmax, achieves 3/4. This fact indicates 
that if the AN z is properly generated and Ez is above a certain 
threshold, applying ideal coherent demodulation at the relay 
does not lead to better performance than the non-informative 
case, and secure data transmission can therefore be ensured by 
the proposed AN design. Eigure 9 shows the SER performance 
for square 16-QAM constellations. Comparing Eigure 9 with 
Eigure 8, similar observations can be made as for the case of 
4-QAM. 

Now, we investigate the ASER performance at the relay 
using the optimal AN design in Section IV. Eigure 10 depicts 
the ASER performance for square 4-QAM and 16-QAM 
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Fig. 9. SER performance of 16-QAM modulation as a function of Ez/N q, 
where Em/No = 10 dB, h = 1, and g = 1. 



Fig. 10. ASER performance of 4-QAM and 16-QAM modulations as a 
function of Ez/Nq, where Em/No = 10 dB, = 1, and Ug = 1. 

signals. In this figure, ASERmax is computed by numerically 
maximizing ASER(a;i,X 2 ) in (26) with respect to xi and X 2 - 
Before comparing with other power distributions, we note that 
from numerical results, the optimal power allocation for the 
considered scenarios in Figure 10 is ^ with probability 
p = 1. In other words, different from the AN design based 
on the instantaneous CSI in Section III, a constant power 
usage is optimal to design the AN based on the statistical CSI. 
Intuitively, due to the averaging effect over several channel 
realizations, the ASER performance is not sensitive to the 
different instantaneous power allocations. For comparison, the 
ASER curves for various power distributions are plotted, i.e., 
the uniform and the exponential power distribution. For each 
curve, the numerical results are obtained from computing 
ASER(Y^)/(a;)dx, where 

fix) = for 0 < x < 2P (27) 



f indB 


Fig. 11. ASER performance of 16-QAM modulation as a function of Ez/Nq, 
where Em/No = 10 dB, (Jh = 1, and <7g = 1. 

and 

/(a;) = iexp(^-p^ , fora;>0 (28) 

are the PDFs of the AN power with uniform and exponential 
distributions, respectively. Note that there exists numerous 
AN designs with exponentially distributed power, such as 
the Gaussian distributed AN. Figure 10 shows that our pro¬ 
posed scheme with the SER-maximizing power distribution 
yields the largest ASER at the relay. Given sufficiently high 
ANR, all the ASER curves achieve the non-informative error 
performance, which confirms Proposition 5. However, the 
performance difference between different AN designs is not 
significant. This can be explained by the fact that the ASER 
expression in (25) does not consider the phase of the AN. 
Therefore, the degree of freedom of the phase design for the 
AN is not utilized. 

Furthermore, it is interesting to compare the SER perfor¬ 
mance of the optimal AN designs in Sections III and IV. 
Figure 11 compares the ASER performance of various AN 
designs. In the figure, “Instantaneous CSI Based Sim.” curve 
is plotted by computing f ft, g [SERmax] using the Monte Carlo 
method, and the curve “Statistical CSI Based Ana.” is plotted 
using the same method as in Figure 10. The curve “Gaussian 
Sim.” is plotted by numerically computing fft g[SER(z)], with 
Gaussian distributed AN. From Figure 11, we observe that the 
instantaneous CSI based AN design yields a higher ASER than 
the statistical CSI based AN design, which is consistent with 
our intuition. However, the performance difference between 
these two designs is not quite significant. Without optimal 
phase and power designs, the Gaussian distribution performs 
worse than both CSI based AN designs. 

VI. Conclusions 

In this paper, we have investigated physical layer secrecy 
for a two-hop single-antenna relay channel, where one source 
aims to transmit to one destination assisted by one untrusted 
relay. In order to achieve physical layer secrecy, a novel AN 
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scheme was proposed to enhance the secrecy of wireless com¬ 
munications from the source to the destination. The objective 
was to design AN symbols generated by the destination node 
that degrades the error probability performance at the untrusted 
relay, in particular, that maximizes its achieved SER perfor¬ 
mance. For the case where perfect instantaneous CSI of the 
source-to-relay and relay-to-destination links is available at the 
destination, we have derived exact analytical SER expressions 
for the relay and studied the optimal design of the AN signal 
to maximize the corresponding SER. In this work, the source 
adopts a squared M-QAM modulation scheme as in current 
modern cellular standards. It was shown that the Gaussian 
distribution, which is frequently used in the context of AN 
[8,50,51], is not optimal in general. The optimal AN for the 
considered relay channel was found to be QAM or rotated 
QAM phase selection. Moreover, compared with the Gaussian 
AN, our optimal AN distribution can yield remarkably higher 
SER at the relay. For the case where the AN design is 
based on the long-term CSI of the source-to-relay and relay- 
to-destination links, we presented the corresponding ASER 
expression at the relay, and the optimal AN distribution to 
maximize the ASER performance was determined accordingly. 
Interestingly, it was shown that the phase of the AN does 
not affect the ASER performance at the relay. Regarding the 
design of the power distribution of the AN, the obtained 
power distribution was shown to deliver improved ASER 
performance than various other power distributions, such as 
the uniform and exponential distributions. Finally, we mention 
that there are many directions for further extensions of this 
work. For example, when perfect CSI of the source-to-relay 
and relay-to-destination links is not available, opportunistic 
scheduling with low-rate CSI feedback might be applied 
[52]. Other interesting directions include the AN design for 
imperfect CSI knowledge as well as for the multiple-antenna 
relay channel and two-way communication networks. 

Appendix: Proof of Theorem I 

In the following, we base our proof on the concept of Linear 
Programming and the Karush-Kuhn-Tucker (KKT) optimality 
conditions [53, p. 243]. 

Before proceeding, note that for a discrete realization xi 
with probability p, one can express f{x) at xi as 

f{x)=pd{x-xi). (29) 

If a PDF f{x), which solves the problem (20), has the form of 
f{x)=S{x — xo), we can easily obtain xq = P, as SER (x) 
in (19) is a monotonically increasing function. This case is 
trivially contained in Theorem I. In the following, we will 
assume that /(x) > 0 for at least two different values of x. 
In this case, at the optimum, the constraint (20b) must be 
met with equality as otherwise we can further increase the 
objective function ( 20 a) by increasing f(xi) and decreasing 
f{x 2 ) for some xi > X 2 without violating any constraint, 
which contradicts to the optimality assumption. 

If the optimal distribution f{x) corresponds to that of a 
discrete random variable with n realizations, f{x) can be 


expressed as 

n 

f{x) ='^p^5{x-Xi) (30) 

i=l 

where Pi > 0- Inserting (30) into (20a)- 

( 20 c), we can obtain the following system of linear equations 


aipi + a 2 P 2 H-f anPn = b (31a) 

PlPl + P 2 P 2 H-f PnPn = P (31b) 

Pl+P2-\ - VPn = l (31c) 

where = SER(yx 7 ), A = Xj, and 0 < < 1 with 

i = 1,... ,n. The achieved maximum value of the objective 
function ( 20 a) is denoted as b. 

If the optimal PDF /(x) contains at least one interval 
[xa,Xf,] with /(x) > 0 for any x G [xo,Xf,], we can divide 
[xa,Xf,] into m non-overlapping sub-intervals {xm-i,Xm) 
with Xa = xo < xi < • • • < Xm = Xb. Then, following 
the first mean value theorem [54, Theorem 12.111], we can 
find some Xi and Xi with Xi-i < Xi < Xi and Xi-i < Xi < Xi 
such that 

Pi= [ fix) dx (32) 

J Xi_i 

m 

SER {\/x) f (x) dx = SER Pi (33) 

and 

m 

xf{x) dx = ^ XiPi. (34) 

Denoting SER (v^) and Xi as at and /3j in this case and 
repeating the approach for all intervals on which /(x) > 0, 
we can again obtain a system of linear equations in the form 
of (31). Note that Xi —>■ Xi-i as the number of sub-intervals 
m increases. In this case we have Xi ~ Xi ~ Xi_i ~ Xi-i. 
Further, as SER (x) is a monotonically increasing function, 
we can order the coefficients such that ai < a 2 < ■ ■ ■ < o-n 
and Pi < 132 <■■■ < Pn- Then, as 0 < pi < 1 for all 
i = 1, • • • , n, (31b) and (31c) indicate that Pn > P. 

As /(x) is the optimal distribution and b is the maximum 


of the objective function (20a), pi with i = 1, 2,..., n in (31) 
must be a solution of the following problem 

max ctiPi -I- a 2 y 2 H - f a„p„ (35a) 

S.t. Piyi + P2y2 H-h PnVn = P (35b) 

2/1 + P 2 H-f Pn = 1 (35c) 

0 < Pi < 1 for / = 1, 2,..., n (35d) 




where and Pi with i = 1 , 2 ,..., n are identical to the cor¬ 
responding parameters in (31). Following the KKT conditions 
[53, p. 243], at the maximum of the problem (35), we must 
have 


Ao = 


ai — Uk 

Pi - Pk 


(36) 


for any i ^ k with 0 < Pi < 1 and 0 < p^ < 1 , where Aq is 
a constant (Lagrange multiplier). Without loss of generality, 
let us denote k' as the first index such that 0 < Pj < 1 , i.e.. 
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k' < i for any 0 < pk',Pi < 1- Similarly, we can denote 
k* as the last index such that 0 < pi < 1, i.e., k* > i for 
any 0 < pk*,Pi < 1- Clearly, we have Pk' < P < Pk*, as 
otherwise pi with i = cannot meet the constraints 

(35b)-(35d). Inserting (36) in (35) and eliminating all ai and 
Pi for i ^ k', the optimum value of b can be expressed as 

b = \qP + ak’ — XoPk' ■ (37) 

Moreover, let 

yk' = -^—and y, = 0 for i p kPk* 

Pk- — Pk' Pk* — Pk' 

(38) 

and consider the following PDF f{x) 

Inserting (39) into (20), we obtain the same maximum value b 
for the objective function (20a) as given in (37) while meeting 
all the constraints. Therefore, there exists a PDF f{x) solving 
the problem (20) with f{x) > 0 for at most two values of 

X. ■ 
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